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Abstract. We examine the exact perturbations that arise from the q-average 
formalism that was applied in the preceding article (part I) to Lcmaitre-Tolman- 
Bondi (LTB) models. By introducing an initial value paramctrization, we show 
that all LTB scalars that take a FLRW "look alike" form (frequently used in the 
literature dealing with LTB models) follow as q— averages of covariant scalars that 
are common to FLRW models. These q— scalars determine for every averaging 
domain a unique FLRW background state through Darmois matching conditions 
at the domain boundary, though the definition of this background docs not require 
an actual matching with a FLRW region (Swiss cheese type models). Local 
perturbations describe the deviation from the FLRW background state through 
the local gradients of covariant scalars at the boundary of every comoving domain, 
while non-local perturbations do so in terms of the intuitive notion of a "contrast" 
of local scalars with respect to FLRW reference values that emerge from q- 
averages assigned to the whole domain or the whole time slice in the asymptotic 
limit. We derive fluid flow evolution equations that completely determine the 
dynamics of the models in terms of the q-scalars and both types of perturbations. 
A rigorous formalism of exact spherical non-linear perturbations is defined over 
the FLRW background state associated to the q-scalars, recovering the standard 
results of linear perturbation theory in the appropriate limit. We examine the 
notion of the amplitude and illustrate the differences between local vs non-local 
perturbations by qualitative diagrams and through an example of a cosmic density 
void that follows from the numeric solution of the evolution equations. 
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1. Introduction. 

in the preceding article (part I) we introduced for the study of LTB models pQ a 
formalism based on a new set of scalar variables (the q-scalars) that follow from 
applying a weighed proper volume average (q-average) to covariant fluid flow LTB 
scalars that are common with FLRW models. As proven in part I, the q-scalars 
are not coordinate ansatze, but covariant scalars related to curvature and kinematic 
invariants, and thus provide an elegant and coordinate independent representation of 
the models that is alternative to the standard variables normally used in the literature 

0[31H|5JH[71[E1[H] ( see reviews in EB EH! ED El])- All proper curvature and 
kinematic tensors characteristic of the models are expressible in terms of irreducible 
algebraic expansions formed with the metric and 4-velocity, whose coefficients are local 
fluctuations of these scalars. Also, all scalar invariant contractions of these tensors 
are quadratic fluctuations of the q-scalars whose q-averages are statistical moments 
(variance and covariance) of the density and Hubble scalar expansion. 

As shown in part I (see summary in section 2), the q-scalars can be, either 
functionals defined on arbitrary fixed domains, or functions ( "q-functions" ) when 
considering the pointwise dependence of the average on the varying boundary of a 
domain. By comparing q-scalars with the non-averaged covariant "local" scalars 
we obtained fluctuations and perturbations (see section 2), which are exact, not 
approximated, quantities. The fluctuations and perturbations can be "local" when the 
comparison is with q-functions in a pointwise manner, or "non-local" if comparing 
local non-averaged values with the q-average assigned to a whole domain. 

The q-functions and their corresponding local perturbations have been applied 
successfully to examine various aspects and properties of LTB models: to construct 
an initial value formulation [13] . to examine inhomogeneous dark energy sources 
(quintessence and the Chaplygin gas) [T^l [T5] . to apply a dynamical systems approach 
[HH H7) . to examine their radial asymptotics [TS], the evolution of radial profiles 
of covariant scalars and void formation 19 , to probe the application of Buchert's 
averaging formalism to LTB models [20l [2lJ [22] and to study the dynamics of non- 
spherical Szekeres models [53] . In the present article we extend and enhance previous 
work by considering also non-local perturbations and by discussing various properties 
of all perturbations not examined previously (their extension, amplitude, their use in 
Swiss cheese models and asymptotic properties). 

As a continuation of part I, we examine the q-scalars and their perturbations 
(local and non-local) in the framework of an initial value parametrization that is 
introduced in section 3, so that all relevant quantities can be scaled with respect 
to their values at an arbitrary fiducial (or "initial") time slice. This initial value 
parametrization emphasizes the role of q-scalars as LTB scalars that behave as 
"effective" FLRW scalars, as they (i) satisfy FLRW time dynamics, (ii) mimic FLRW 
expressions that are widely used in the literature (for example, in many of the void 
models Pl[5Sl[2i[571[5Sl|23[3J[ni[^[33IMI^),and (iii) identify for each domain 
a unique FLRW background state through Darmois matching conditions (though an 
actual matching with a FLRW region is optional, not mandatory). 

Considering that (as proven in part I) local perturbations convey the deviation 
from FLRW geometry through the ratios of Weyl to Ricci curvature invariants and 
anisotropic (shear) to isotropic expansion (see equations I(42a)-I(42b) and 1(43) [f] 

X We will frequently use equations derived or presented in part I. We will refer to these equations by 
the notation "I(X)", where "X" corresponds by the equation number in part I 
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), and bearing in mind that these perturbations and their associated q-functions 
completely describe all proper tensors and scalar contractions, it is natural to expect 
that these q-scalars and their perturbations should also yield a complete and self- 
consistent system of evolution equations that fully determine the dynamics of the 
models [HI [13 H3 H3 [TSJ, \T§\ . In section 4 we derive these evolution equations for 
local and non-local perturbations. 

The evolution equations for the q-scalars and their perturbations (local and non- 
local) have the structure of evolution equations for exact spherical perturbations over 
the FLRW background (now described in terms of the q-scalars through Darmois 
matching conditions). Since q-scalars are covariant LTB objects that satisfy FLRW 
scaling laws and time evolution equations, while the perturbations (local and non- 
local) effectively convey the deviation from FLRW behavior, it is natural that a 
rigorous formalism of exact spherical perturbations on a FLRW background can 
emerge from this set of variables in which the q-scalars common to FLRW are 
"zero order" variables and the remaining scalars (including local scalars) are "hrst 
order" quantities. As shown in jT^l [TB] , such formalism arises for the case of local 
perturbations. We extend in section 5 this result to non-local perturbations. 

An important distinctive characteristic of the fluctuations and perturbations 
(local or non-local) is their extension along the radial range: they can be either 
"confined" (i.e. localized in a given bounded comoving domain, with or without a 
matching with a FLRW region as in a Swiss cheese model) , or "asymptotic" (when the 
domain becomes the whole time slice). We discuss the difference between confined and 
asymptotic perturbations in section 6, showing (in particular) that local perturbations 
can always be treated asymptotically, whereas asymptotic perturbations can be non- 
local only for LTB models that converge to a FLRW model in the asymptotic radial 
range [TH]. In this latter case, the perturbations measure the "contrast" of local 
scalars A with respect to a global reference value given by the asymptotic limit of 
the q-average functional, which coincides with the equivalent scalar A of the FLRW 
asymptotic state. This type of asymptotic perturbations is often used in the literature 
when considering perturbations in the context of LTB models (see examples and 
reviews of these "contrast" perturbations in [TTJ Q2], see also the linear regime in 
the Appendix of [3l)]L 

Since local and non-local perturbations (whether confined or asymptotic) are 
different objects, we illustrate this difference in section 7 by showing how they provide 
a different measure of the deviation from the FLRW background: local perturbations 
describe this deviation through the local magnitude of radial gradients of covariant 
scalars, whereas non-local perturbations describe it through the familiar and intuitive 
notion of a "contrast" between local values of non-averaged scalars and reference 
average values assigned to a whole domain (or to a whole slice in the asymptotic 
limit). As a consequence, non-local perturbations are more intuitive, as the sign of 
their amplitudes corresponds to the familiar positive/negative sign that we associate 
to over/under densities. This sign is the opposite for local perturbations: over/under 
densities are negative/positive. 

In spite of their differences, we show in section 8 that in the linear limit both 
perturbations (local and non-local) yield the familiar density perturbation equation 
of linear theory in the comoving gauge. In section 9 we use a numerical solution of the 
evolution equations derived in section 4 to present an example of a cosmological void 
configuration that converges to an Einstein de Sitter FLRW model in the asymptotic 
radial range. Besides being useful to appreciate the difference between local and 
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non-local asymptotic perturbations and their connection to the radial profiles of LTB 
scalars, this numerical example shows the potential utility of the q-scalars and their 
evolution equations for LTB model construction. We provide a summary and final 
discussion in section 10, while in Appendix A we prove that shell crossing singularities 
necessarily occur in Swiss cheese models matching a hyperbolic LTB region with an 
Einstein de Sitter background. 



2. The q average, q scalars and their perturbations. 



LTB dust models are characterized by the metric 1(1) and the field equations I(2a) 
I(2b), which we repeat below for convenience: [f] 

J? 2 
1 + 2E 



ds 2 = -dt 2 + — dr 2 + R 2 (dtf 2 + sin 2 dp 2 ) , (1) 



R 2 = ^ + 2E, (2a) 

2M' = %TtpR 2 R\ (2b) 

where R = R(t,r), R = u a V a R = dR/dt, R' = dR/dr, E = E(r),M = M(r) and 
p = pit, r) is the rest mass energy density (we have set G = c = 1 and r has length 
units). The basic covariant fluid flow scalars of the models 1(7): 

6 

p, % = - (Hubble expansion), /C = — — (spatial curvature), (3a) 
3 6 

£ (eigenvalue of the shear tensor), £ (eigenvalues of the electric Weyl tensor) (36) 



2M 



where the rest mass density p is given by (26), the expansion scalar is 6 = V a u a and 



TZ^ is the Ricci scalar of the hypersurfaces 3 T[£] orthogonal to u a (the time slices) 



The local scalars pd)-(36) can be computed from the metric functions by means of 



(26) and 1(3)— 1(6), and their 1+3 evolution equations and constraints are given by the 



system I(8a)-I(8d) and I(9)-I(10). 



2.1. The q-scalars. 

The LTB scalars that are common to FLRW spacetimes are ( |3a[ ): p, H, K., and as we 
showed in part I, their q-averages (defined by 1(13)) in an arbitrary fixed spherical 
comoving domain 2?[r&] are given by the functionals 1(14)— 1(16) : 

f(p) 9 N = §, (n) q [r b } = ^, </c> g M = -||, (4) 

that satisfy the constraint 1(17): 

{U) 2 q [n] = y<p) 9 N - {K) q [r b }, (5) 

where the subindex /, indicates evaluation at r = r&. The functionals above assign the 
real numbers in the right hand sides of Q for the whole domain T>[ri)]. However, if we 

§ This section provides a quick summary of results of the preceding article (part I) that will be 
needed in the present article. For more detail and explanation the reader is requested to consult part 
I. 
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consider the q-average definition 1(13) to construct real valued functions depending 
on a varying domain boundary, we obtain q-functions: p q , T-L q , K, q that comply with 



4tt M R „ 2E 

8tt 



which is formally identical to Q and ([5]) but hold in a point-wise manner for every r 
(see [213 [211 122] and part I for more detail on the difference between the functional 
(A)g[rh] and functions A q (r)). Notice that the q-scalars (either as functions or as 
functional are not coordinate '"ansatze, but fully covariant objects since M and R 
are invariant scalars in spherically symmetric spacetimes [37] and E is related to them 



via (2a) (see part I for further discussion on this issue). 

Since quantities that are functions of q-scalars are themselves q-scalars (see 
Appendix B of part I and Appendix B of (23j), then we can define the q-scalar 1(25) 
= U(pq,H q ) or fl q = U{p q ,K, q ) either as a functional or as a function by: 



which is formally identical to the FLRW Omega factor. It is straightforward to show 
that the q-scalars p q , Ji qi IC q , Q q (whether evaluated as q-functions or as functionals 
in fixed domains 2?[r(,]) satisfy the FLRW evolution laws I(27a)-I(27b). fjjl which 
evidently single out the q-scalars as LTB scalars that behave as FLRW scalars (in the 
sense that they comply with FLRW time dynamics). 

2.2. Local perturbations. 

If A and A q = (A) q are both evaluated as real valued functions on the same arbitrary 
value r that denotes a varying boundary of concentric domains T>[r] for r > 0, then 
a local perturbation follows by the pointwise evaluation comparison at each r of the 
ratio 1(29): 

r(A) M _ Mr)-(A) q [r] _ A(r)-A q (r) 

[) (A) q [r] " A q (r) ' W 

which comply (from 1(23) and 1(24)) with 1(30) that relates the 5^ with radial 
gradients of A q and A (also valid for the (A) q ): 

that leads, using 1(23), 1(24) and I(B3), to the following linear algebraic relations 
among the 5^: 

25™ =n q 5 {p) + [l-tt q ]5 {k \ (10) 

<y(") = got - 2sm = (i - n,) (>> - s^A , (11) 

where Q q is given by til and 8^ above is consistent with £1 = fl q (l + 6^) in 1(26). 



|| For the functionals (A) 9 [r(,] the derivatives involved are (A)'q [rj] , which can be evaluated (at r = ri) 
either directly from (HJ, or with the commutation rule 1(22) and the forms of the local (non-averaged) 
scalars in \2b\ , 1(3), 1(4) and (JTl. If using 1(22) for computing ("H}^;,] and (f2)' q [ri,] we also need to 
use the identities 1(36) and 1(37) that are proved in Appendix C of part I. 
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2.3. Non-local fluctuations and perturbations. 

As opposed to local perturbations in which A and A q = (A) q evaluate at the same r, 
we can define for every fixed domain T>[r b ] non-local perturbations 

^Hr,r b ) = A{r) { -^ [rb \ 0<r<r b , (A) q [r b }^0, (12) 

that compare local values A(r) inside the domain with the q-average (functional) of A, 
which is a non-local quantity assigned to the whole domain (notice that at every 3 T[t] 
the value (A) q [r b ] is effectively a constant for all r < r b and a function of t for varying 
3 T[t}). Evidently, the 5j$(r, r b ) do not comply with |9| and the properties that follow 
thereof (notice that d/dr[S^ (r, r b )] = A' / {A) q [r b ]). As shown in part I, the non-local 
fluctuations D NL (A) that give rise to non-local perturbations (D NL (A) = 8i^}{A) q ) 
are effectively statistical fluctuations. 



3. The q scalars define a FLRW "background state". 



3.1. The q-scalars as LTB objects that look like FLRW expressions. 

"FLRW look alike" expressions are often introduced in various applications of 
LTB models, specially in a lot of recent articles looking at LTB void models 
[241 [251 EH EH EH E9J EQl EH E21 [331 EH E5]. These expressions are introduced in 
these references as "convenient" ansatzes, without any justification other than their 
"FLRW look alike" forms, thus ignoring the fact that they can be defined rigorously 
as q-scalars that emerge from the weighed average 1(13), and thus are fully covariant 
quantities related to curvature and kinematic invariants (see section 6 of part I) . The 
"FLRW look alike" expressions follow readily by parametrizing the metric functions 
in (JlJ in terms of their fiducial values at an arbitrary time slice t = t . Considering 
the coordinate choice R = r [^] , we can transform (JlJ into the FLRW "look alike" 
metric: 



ds 2 = -dt 2 
R 

Rq 



Fdr 



1 - K. q0 r 2 



+ r 2 (dt? 2 + sin 2 M<p 2 ) 



R 



= H 



Y = R '/ R 



Rq/Rq 



rR 1 
~R 



= 1 



(13a) 



(136) 



where the relation between the scale factor a and % 



follows from (66 1 (see section 
(fl~36| the q-scalars p q , 



7 of part I). Under the parametrization (13a 

(6a)-(66) and ^ (and their functional equivalents) take the following "FLRW look 
alike" forms that often appear in the literature: 



Pq 



PqO 



n q — H q 



n 



a 

1 - Q 



q o 



n q = 



n 



qO 



ClqO + (1 - Mqo) a 



(14a) 
(146) 



where the subindcx q indicates evaluation at t = to. 



1 This coordinate choice is not appropriate for LTB models whose time slices have spherical § 3 
topology or lack symmetry centers. For such models R' (and thus R') are no longer monotonical on 
r: they must change sign at a fixed r (a turning value) at every time slice. This turning value is a 
common zero with the gradients of all scalars |19) . 
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3.2. The FLRW "background state" through Darmois matching conditions. 



For every fluid flow FLRW scalar A = p, H, K,, fl (we denote henceforth FLRW scalars 
by a tilde) and its "FLRW equivalent" LTB q-scalar A q = p q , H q , JC q , Q q , the value 
Agb = A q (t,rb) identifies, for each domain D[r&] of an LTB model, a specific FLRW 
dust model that could be smoothly matched (under Darmois matching conditions) at 
an arbitrary finite comoving radius r — rb that also marks the boundary of 2?[r&]. 
Consider a dust FLRW universe with metric 

d;J 2 (di? 2 + sin 2 tfd^ 2 ) 



As 1 = -&t z + a\t) 



(15) 



l-k Q q 2 r* 

where k = 0, ±1 , £q is an arbitrary length scale and a(t) is the dimensionless FLRW 
scale factor. Darmois conditions (necessary and sufficient) for the smoothness of the 
matching of ( 15 ) with an LTB model along a comoving boundary r — rj, are given by 

mm 



Pqb = {p)q[n} = P{t) 



PO 



K-qb = (IQq[rb\ = lC(t) 



(16a) 



n qb = (H) q [n] =u(jt) = -, ci qb = (n) q [r b ] = n(t) = 



8irp 

W 2 



Rb = a(t) r or a& = a, 

where the subindex denotes evaluation at a fiducial hypersurface t 
holds. We note that the continuity of the q-scalars under the matching conditions 



(166) 
(16c) 

to and a(t ) = 1 



( 16a 



(14a 



and (16b) is strikingly evident if we use the parametrization ( 13a)— ( 136 1 and 
-(146) with Aq = A q t>o — A q (to,ri,) holding for A — p, %. IC, ft. However, from 



(91) and (3a)-(3b), it is evident that Ab ^ A q b and A' q , A' q ^ hold in general, and 

"36| and the gradients A' and A' q do not comply with the 
q). 



thus the local scalars (3a 



matching conditions ( 16a 



It is important to remark that the continuity of the A q under Darmois matching 
conditions is simply a formal rigorous procedure to identify for every T>[rb] of an LTB 
model a particular FLRW dust model that can be defined as a reference "background 
state". We use the term "state" to emphasize that this identification does not 
force us to consider an actual matching with a FLRW region, which would yield a 
"Swiss cheese" configuration in which the background state becomes also an actual 
background spacetime. Likewise, the A q also allow us to define a "FLRW equivalent" 



region to every domain 2?[r(,], as they provide through (16a)-(16c) the values of the 



FLRW scalars A if the whole domain was replaced by an equivalent spherical comoving 
section of a FLRW spacetime (without mass or volume compensation [55]V 

Therefore, having determined the FLRW background state, the discontinuity of 
the local scalars and the gradients A' and A' q is not problematic if we do not wish to 
construct an actual Swiss cheese model through a smooth matching at the domain's 
boundary r = rb- In this latter case we can avoid discontinuities by demanding (besides 
(16a)-(16c)) also the continuity of local scalars at r = r& through the following extra 



supplementary condition: 

A' q b = 0, 



=> A h = A 



qb 



for A = p, K, H, fi, 



(17) 



which forces A to coincide with A q at r = rb, and thus explains the appearance (see 
panels (b) and (d) of figures 1 and 2) of "humps" (if A' > 0) or "bags" (if A' < 0) in 
the radial profiles of local scalars (this has been noted in the local density profiles in 
Swiss cheese models in LTB void models [301 EI] , see reviews in chapter 5.3.5 

of QU and in [35]). 
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4. Evolution equations. 

While the q-scalars A q or (A) 5 [r&] satisfy FLRW evolution laws, such as I(27a)-I(27b), 
these scalars are not fully determined by these FLRW equations because, unlike their 
equivalent FLRW scalars A, they either depend directly on r or r&. The missing 
dynamical information is provided by the evolution equations for the 5^ and the 

§ (A) 



4-1- Local perturbations. 



As shown in part I, LTB tensors and scalar invariants, as well as the local covariant 



scalars (|3a 
P = Pq 



1 + 5^ 



(36) are expressible in terms of q-scalars and their local perturbations 
H = H a [l + S^-A , K = K a \l + <5 (fe) 



£ 



4tt 



Pq 



(18a) 
(186) 



the evolution equations for the variables p q , H q , S^ p \ S^' should yield a self- 
consistent and complete set of evolution equations for the LTB models. This system 



follows readily by inserting (18a) and (186) into the "1+3" system I(8a)-I(8d) and its 



constraints 1(9)— 1(10) . The result is the evolution equations 

Pq 



ti q 

S(p) 



— 3p q H q , 
n,2 47T 

— n q - 



'Pq, 



3(1 



m = (1 + 3<^>) H q 5W + %^(5M - 

Sri a 



plus the algebraic constraints 
8tt 



■t/2 

H q — 



-Pq - £ 



28M = n q 6<p) + (1 - n fl ) 6^, = 6^ - 2<J< W \ 



(19a) 
(196) 
(19c) 

(wd) 

(20) 



that exactly coincide with the general relations (66), (10) and (11), hence they hold 



at all t (i.e. they propagate in time). The following points are worth remarking: 

• The constraints 1(9) of the 1+3 system I(8a)-I(8d) are satisfied trivially from Q 
applied to <5 (p) and 5^. 



The first two constraints in (20 ) follow from substituting (186) into 1(10) and using 
(2a) and (6a), while the third one is obtained by differentiating Q with respect 
to r and applying Once (19a|-(19dl is solved these constrains and ([7| allow 



us to compute the remaining q-scalars ad their perturbations: fl qi JC q , <5™, S^ 1 ' 



The fact that the constraints (20) of the system (19a)-(19d) are algebraic implies 



a great simplification of the numeric treatment of these fluid flow evolution 
equations, as they can be effectively integrated as a system of autonomous ODE's 
in which the initial conditions are restricted by the algebraic constraints, and thus 
it is far easier to handle than the 1+3 system I(8a)-I(8d) in which the constraints 
1(9) are partial differential equations on r that must be solved before the time 
integration of the equations. The example of the void model given in section 9 



has been obtained from the numerical integration of (19a)-(19d) 
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Considering the relation between the A q and in the constraints (20), each of the 
following combination of four scalars: 

any two of A q = {p q , H q , K. q , fl q }, plus any two of 5 (A) = {5 {p) , 6^, <5 (0) }, 

provides a full covariant scalar representation for the models, since the remaining 
pairs of A q and 5^ can be obtained from these algebraic constraints. Each of these 
representations yields a self-consistent and complete set of evolution equations that is 
alternative (and equivalent) to the analytic solutions of ( |2a[ ) (see Appendix A of part 
I) and to the numerical integration of the "1+3" system I(8a)-I(8d), and thus, they 
completely determine the dynamics of the models. 
The 



evolution equations (19a 



-(19d) correspond to the representation 
{p q , H q , <5 (w) }, which is useful to compare with the spherical collapse model 
and perturbative scenarios of structure formation that consider the density and Hub- 
ble velocity as dynamical variables. However, a more appropriate representation 
for cosmological applications (for example void models) is furnished by the scalars 
{fl q , H q , S^ n \ 6™'}, leading to the following evolution equations: 



H q = -(i + -n q )H 2 q , 



-n q (i-n q )H q 



m = - 



(l + 36W)6W + ln q (5M+5M) 



H 



<j< n > = - [(1 + 36^)6™ - n q (s (H *> + S^) 



H 



(21a) 
(216) 
(21c) 

(21 d) 



though, as opposed to the local scalars p, H, JC, the local scalar ft defined by 1(26) 
lacks a simple direct physical interpretation, besides being a generalization of the 
FLRW scalar f2. 



4-2. Non-local perturbations. 



If we consider non-local perturbations (12 1 in an arbitrary fixed domain 25 [r&], then 
only the local scalars (3a) that are common to FLRW can be expr essed in terms of 



the variables {(A) q [rb\, Si 



p = (p) q [n] 



1 



?(/>) 



} in a similar manner as in (18a) — ( IS 6 1 
1 



H = (H) q [r b ] 



K. = (K)q[r b ] 



1 + 6 {!C) 

1 T U NL 



(22) 



Pq( r ) ~ (P)q[n}5- 



(p) 



(23) 



The remaining local scalars, X and £ cannot be expressed as (36) purely in terms of 
non-local perturbations: 

E(r) = H q (r) - (H) q [r b ](l + 6™), S(r) = ^ 

and as a consequence, the variables {(^4) g [r;,], '}, A = p, H, JC do not provide a 
complete scalar representation of the dynamics of LTB models, which is not surprising 
because for any fixed domain the (A) q [r b ] depend only on t and constitute boundary 
conditions for the A q at r = r b . From (12) and (22), local and non-local perturbations 
are related by 

Mr) 



(A) q [r b ] 



(1 + <5 (A) ), 



(24) 



which shows that they only coincide at the boundary of each T>[r b ] where A q {r b ) = 
(A) q [r b ]. Since their time derivatives are evidently different, it is reasonable to expect 
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that the evolution equations of the 8^ and S^-J will be different. Inserting (I24J) for 



A = p, H into (19c)-(19d) yields these equations: 



Hp) _ 

"NL — 



6™ = (1 



3(l + 8^)(H) q [n}5 { J?\ 
36W)(H) q [r b ] 4l } 



47r (^M ,r(«) 



2(H), [r b ] 1 



<l(H q (r)~(H) q [r b ])6 { *\ 



(25a) 



(256) 



which, in order to render a fully complete system to describe the dynamics of the 
models in an arbitrary fixed domain 2?[rt], needs to be supplemented by the evolution 
equations for (p) q [r b ], (H) q [r b ] and p q , H q : 



(ny q [n] 

(p)' q [n] ■ 



47T 

i q \n\ - y 
-m) q [r b ]{p) q [r 



(p) q [n}, 



Pq\ 



n q (r) = -H 2 q {r) 
l = -m q {r) Pq {r), 



4tt 



(26a) 
(266) 



while the constraints take the form (20) and with the 8^ expressed in terms of the 
by (24) (there is an extra constraint given by (|5j)). It is interesting to remark 



that (25a) is identical to (19c), but (256) differs from ( 19c ) by the terms with H 



which explains the need to add the evolution equations for p q and Ti q . Notice that the 
non-local evolution equations become identical to the local ones (19a)-(19d) at the 
domain boundary in the limit r —> r b where W q (r b ) — (H) q [r b ]. Also the non-local 
evolution equations depend on averages p q (r) = (p) q [r], H q (r) — ( H) q [r ] for inner 



points of T>[r b ] (r < r b ) and are considerably more complicated than (19a)-(19d). 

2661) can be proved easily by comparing 



The self-consistency of the system (25 a 
the mixed derivatives [j^]' obtained from 



[8^'}' = A' (r) I (A) q [r b ] with those that 
follow from the radial derivative of the right hand sides of (25a)-(256|. Evolution 
equations for non-local perturbations on the representation {Q q , Hq, 8 {il \ 8^} can 
also be constructed, but the resulting equations are more cumbersome than (|256| 



and (266). However, it always possible (and easier) to solve the evolution equations 
for local perturbations (either ( |19a |-(( [l9d[ ) or (21a)-((21d)) and then compute the 
non-local perturbations through the relation (24) (this is what was done in section 9). 



4-3. Evolution equations without back-reaction. 

The evolution equations constructed with q-scalars and their perturbations (local 
and non-local) lack the back-reaction correlation terms that appear in the evolution 



equations that follow from Buchert's formalism (notice that (196) is identical to 1(52) 
with Q q = 0) in 1(53)). Also, these equations (in any q-scalar representation) 
form complete and self-consistent systems that can be integrated without further 
assumptions for any given set of consistent initial conditions. On the other hand, 
in order to close and integrate Buchert's evolution equations it is necessary to 
make specific assumptions linking the back-reaction terms with the averaged scalars 
[391 SQl HQ. 



5. A formalism of exact perturbations on a FLRW background. 



It is evident that (as pointed out in previous work [THUS]) the system (19a)-((19d) 
has the structure of evolution equations of spherical dust perturbations (the 8^ A >) on 
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a FLRW background (denned by the A q ). While the non-local perturbations were 
not considered i n the s e refe rence s, th e same resemblance to dust perturbation holds 



for systems like (25a)-(25b) and (26a)-(266| involving (A) g [r&], A q (r) and 8i^}{r,r b ) 
Following the standard methodology [36i WR fl3l I3H 1351 14*51 ITfl I"""*"] . a perturbation 
formalism linking a "lumpy" spacetime (LTB model) and a homogeneous "background 
spacetime" (a FLRW dust model) in a domain T>[r] can be defined by comparing local 
LTB variables with LTB objects that can define an FLRW "background state" of 
"zero order" variables by means of suitable maps (evidently, these objects are the 
q-scalars). However, such maps must also deal with gauge issues involving a specific 
time slicing and coordinates. Since the boundary of every spherical comoving FLRW 
region can be mapped into the boundary of a domain T>[r] of an LTB model by the 



Darmois matching conditions ( 16a )—( 166 1 , and bearing in mind that LTB models and 
dust FLRW universes are both (i) spherically symmetric, (ii) have a geodesic 4-velocity 
and (iii) their full dynamics reduces to scalar modes 36, 46], a perturbation formalism 
associated with the q-scalars in domains T>[r] can be defined rigorously by means of 
maps between FLRW covariant scalars and the q-scalars with all gauge issues resolved. 

5.1. The perturbation maps. 

Let Ai be an LTB model and X(T>[r]) the set of all covariant scalars in an arbitrary 
domain T>[r] of Ai. Let Ai be a dust FLRW model and X(Ai) the set of covariant 
scalars of Ai: 

The local perturbation map. For every T>[r] in Ai there exists a model 
Ai such that A q (t, r) = A(t) holds for every A <E X{M) and A £ X(M) and 
for all t. The following maps 

$ : X(M) -> X(V[r}), A h> = A q (r) € X(V[r]), (27) 

> : X (m) -, Xm) , A » = ^> - ^ . *=gk, 

$(A) A q {A) q 

(28) 

define for every T>[r] a "background state" associated with an FLRW 
cosmology M. and local exact perturbations of scalars A obtained by 
comparing them with the LTB scalars produced by the map $. 

An analogous scalar perturbation formalism for the non local fluctuations can be 



defined along the lines of (27) and (28): 



The non-local perturbation map. Let F(Z>[rb]) be the set of all linear 
functionals in an arbitrary fixed comoving domain V[r b ] C 3 T[t\. For every 
FLRW covariant scalar A G X{M.) and every T>\r^\ the following map 

$ NL : X{M) -> Y(V[r b }), A^(A) q [r b ], (29) 

defines a "background state" associated with a FLRW cosmology M 
but consisting of the functionals 1(13). The definition of the non-local 
perturbation of the scalar A is then the map Sil : X{M) x X(T>[r b ]) — > 
X{V[r Q ]) x Y(V[r ]) such that 

$nl(A) (A) q [rb\ 
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Figure 1. Local confined fluctuations. The fluctuations D(A) = A — A g = 
Aq^' 71 ' (thick vertical lines) follow by pointwise comparison of the scalars A and 
A q along r in an arbitrary time slice 3 7~[t]. Clump and void profiles are displayed 
in panels (a)— (b) and ( c)— (d), the case of a Swiss cheese model matched to a 
FLRW region at r = r is displayed in panels (b) and (d). Notice how in these 
latter cases D(A) = at r = rj, thus forcing A q (r^) = A(r(,), which explains 
the "humps" and "bags" in the radial profile of A. Shell crossing singularities 
necessarily emerge in Swiss cheese models in which the LTB region is hyperbolic 
and the FLRW region is Einstein de Sitter (see Appendix A). 



which provides a pointwise comparison between the real valued functions A 
and their associated functionals (A) g [rf,] assigned to the whole domain 2?[r&]. 

It is important to emphasize that in the definitions above we distinguish between 
the FLRW "background state" and the FLRW "background spacetime" M. As 
commented in section 3, the former is defined by the q-scalars: {p q , H q IC q , il q } 
or {(p)q[ r b\, (H)q[ r b], (fc)q[ r b], (tyq[ r b]}, which satisfy FLRW dynamics and relate to 
the FLRW scalars {p, H, K,, f2} of the background spacetime M. by continuity under 
the Darmois matching conditions at an arbitrary constant r. These q-scalars are (in 
the perturbation maps) the "zero order" variables, with the "first order variables" 
being the local fluid flow scalars A related to the former by the perturbations 5^ A ' or 

S (A) 



Notice that the perturbation formalism defined by (27 1 and (28 1 is covariant 
because the background variables and the perturbations (zero and first order variables) 
are coordinate independent LTB objects (see section 6 of part I). In fact, following the 
Stewart-Walker lemma [12 031 0H 05] , the perturbations 5 (A ^ and 5^ are also gauge 
invariant (even in the usual sense of [50]) as they vanish in the FLRW spacetime M. 



associated with the background state through (27 1 and (28) 
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Figure 2. Non-local confined fluctuations. The panels display (as in figure 
1) clump and void profiles with and without a Swiss cheese matching to a FLRW 
region at r = rj, in an arbitrary time slice 3 T[t]. The non-local fluctuations 

Dnl(j4) = A — (A) q = (A) q &\^ (the arrows) compare local values of A(r) with 
the average (j4)g[fi,] (thick line) assigned to the whole domain ©[»"(,] (shaded 
region). Notice (as in figure 1) that D N l(j4) only vanishes at r = ri, if there is a 
Swiss cheese matching (panels (b) and (d)). 



5.2. LTB models as exact perturbations. 

As opposed to the conventional approach to perturbations, either the traditional 
one with gauge invariant variables [UJ HH [5U] or the covariant formalism of Ellis 



et al @2J HH HI]), the perturbations that emerge from (|27)-(|28|) and ([29]) — ([30]) are 
exact (not approximate) quantities, and thus do not lead to some unknown "near 
FLRW" space-time on the basis of a linearization process (though a rigorous linear 
limit can be defined, see section 8). Instead, the {A q , 5^} or the {(A) 9 [ro], 5^} 
with A = p, %, fi, K express a known class of spacetimes (generic LTB models) as 
exact spherical perturbations on an abstract FLRW background state defined by LTB 
objects: the A q along all domains T>[r] with varying boundary or the (j4) 9 [r{,] at a 
fixed domain. Evidently, these perturbation formalisms are special in the sense that 
they preserve the spherical symmetry and the dust source of the FLRW background. 
However, these formalisms are applicable to any space-time compatible with an LTB 
metric in the comoving frame and having an anisotropic fluid source [15) . and can be 
readily generalized for the non-spherical Szekeres dust models [23] . 

6. Confined perturbations, Swiss cheese models and asymptotic 
perturbations. 



In the definitions (27|-(28l and (29)-(30l we assumed domains T>[r] or T>[ri,] with r 
and rf, finite. As we show below, the resulting perturbations can also be considered 
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in the asymptotic case when rt, — > oo. 

6.1. Perturbations confined in comoving domains. 

As long as r and rj, are finite, we are considering perturbations that are confined 
in arbitrary bounded comoving domains of an LTB model. Depending on whether 
we perform a smooth match with a FLRW spacetime or not we have the following 
possibilities: 

There is no matching with FLRW (see panels (a) and ( c) of figures 1, 2 and 3). 
The perturbation describes the dynamics of the domain T>[r] or 2?[r&] with respect 



to a FLRW background state defined by ( 27 ) or &29h , which is a fictitious reference 



FLRW dust model Ai (the background spacetime) whose scalars A match (via 
Darmois matching conditions) the zero order variables A q (r) or (.A) g [r(,]. Notice 
that: 

• The fictitious reference FLRW dust model Ai is necessarily different for 
different domains. 

• First order quantities like A and the gradients A' and A' q (A = p, H, JC, £1) 
are not continuous at the domain boundaries r orrj. See section 3 and panels 
(a) and ( c) of figures 1, 2 and 3. 

Swiss cheese holes smoothly matched to a FLRW region (see panels (b) and 
(d) of figures 1, 2 and 3). The perturbation describes the dynamics of a comoving 
domain P[r&] with respect to a background state that in this case corresponds 
to the actual (non-fictitious) FLRW background spacetime matched at r = r^ 



(under conditions (16a)-(16c|) and extending for r > r\>. Notice that: 

• The resulting spacetime is a compound Swiss cheese configuration consisting 
of an LTB section (confined in 2? [?!,]) described by zero order variables A q 
and their perturbations 5^ or 4l , an d the FLRW background spacetime 
that extends for r > r&. The FLRW dust model Ai is necessarily different 
for different domains 2?[r;,]. 

• Darmois matching conditions only require the zero order quantities A q and 
(A) g [rh] to be continuous at r&, though it is always possible to demand ( fl7| ) 
(as extra conditions) so that local scalars A and the gradients A' and 
(which are first order quantities related to via (j9|) are also continuous 
at rfc. These extra conditions imply that the fluctuations and perturbations 
themselves vanish at (see section 3 and panels b and d of figures 1, 2 and 
3). 

• Swiss cheese models have been considered in the literature in the context 
of fitting observations without resorting to dark energy [251 1301 E] 
(see review in chapter 5.3.5 of [12] and in [35]). However, a shell crossing 
singularity necessarily emerge in Swiss cheese models in which a hyperbolic 
LTB region is matched to an Einstein de Sitter background at finite r& (see 
proof in Appendix A). 

6.2. Asymptotic perturbations. 



The application of the perturbation maps (27|-(28) and (29H30) to the asymptotic 



limit r — > oo or r^ — > oo depends on the convergence of LTB scalars and their 
associated q-scalars along radial rays {i.e. curves x a = £<5", where i = J g rr dr is 
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the proper radial length), which are spacelike geodesies of the metric ([!]) and of the 
slices 3 T[t] with metric h ab . As shown in |18) . £(r) is a positive monotonic function in 
regular LTB models, and thus the asymptotic limit r — > oo corresponds to £ — > oo. 

The radial asymptotic behavior of LTB models is determined by the following 
asymptotic limits for the scalars A = p,H, /C, fl: p~] 

A^ (t) ^ 0, Asymptotically FLRW 
lim A(r) = lim A q (r) = lim (A) q [r b ] = I (31) 



n- 



0, Asymptotically Minkowski. 



These limits correspond to q-averages for increasingly large domains up to the 
situation in which T>[r] or 2?[n,] become the whole time slice 7~[t], and thus, we can 
speak of the q-average of a whole LTB model (instead of the q-average of confined 
domains of an LTB model). We look at the perturbations for asymptotically FLRW 
and Minkowski models separately below. 



Asymptotically FLRW models (see figure 4). As a consequence of (28) and 



(31) local perturbations vanish asymptotically in these models: 

lim^(r) = 0. (32) 



However, (30) and (31) lead for non-local perturbations to the following non- 



trivial nonzero asymptotic limit: 

lim 4?(r, r>) = = ^fe^. r finite, (33) 

r i>- > °° 4x>(*) A(t) 

which keeps their interpretation of a "contrast" between local values A(r) and the 
domain average, with the difference that the average now corresponds to a domain 
that encompasses the whole slice 3 T[i] and coincides with the corresponding 
scalars A^ (t) — A(t) of an asymptotic FLRW background state defined by 



(29 1 . As a consequence, we can rigorously state that the q-average of every LTB 
model of this class is exactly the FLRW background spacetime AL Asymptotic 
fluctuations and perturbations for a model converging to FLRW are illustrated 
schematically by figure 4 and for the numerical example of section 9 by figures 6, 
8, 9, 11 and 12b. 



Asymptotically Minkowski models. By looking at (30) it is evident that 
non-local perturbations 5^ cannot be defined in the asymptotic limit for these 
LTB models, since in general A(r) ^ for finite r and thus the ratio A/(A) q [rb] 
diverges as [?"(,] — > in the limit r/, — > oo. However, local perturbations 5^ 
defined by p8] ) are not affected by the radial convergence to a Minkowski vacuum 
because, as shown in (TS], both A — > and A q — > hold as r — > oo but their ratio 
A/A q is finite in this limit, and thus the 5^ in (28) are well defined. The only 



difference with asymptotically FLRW models is that the 6^ ' tend (in general) 
to nonzero constant values as r — > oo |18) . This applies also to the perturbations 
of the spatial curvature 5^ in elliptic and hyperbolic models radially converging 
asymptotically to a spatially flat Einstein de Sitter model (see the profile of the 
perturbation S^' in figure 12b). 

+ LTB models can only converge in the asymptotic radial regime to FLRW dust models with zero or 
negative spatial curvature, not positive. Models convergent to Minkowski can converge to a section 
of Minkowski parametrized by Milne coordinates or by non-standard curvilinear coordinates. In the 
Milne case p q , £l q — > but H q , K q tend to nonzero values (see 1181 ). 
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7. Comparison between local and non local perturbations. 

Since the perturbations and the 6^ are different objects, they provide a different 
measure of the inhomogeneity or "deviation" of an LTB model from a FLRW Universe: 

• The 5^ measure this deviation through local comparison of A(r) with the 
characteristic FLRW scalar A q (r) = [A) q [r] at every r, and is related (via 

to radial gradients of A q and A and (via I(42a)-I(42b) and 1(43)) to the ratio of 
Weyl to Ricci curvature and anisotropic to isotropic expansion. See figures 1 and 
3. 

• The 

<5ni? measure the deviation from FLRW though the notion of the "contrast" 
between local values of A[r) inside an arbitrary but fixed comoving domain 
(0 < r < r&) and a value (j4) g [r&] that can be associated to the equivalent FLRW 
scalar A(t) = (A) q [ri,] that characterizes the whole domain. In the asymptotic 
limit rf, — > oo so that A(t) = A^ (t), these perturbations provide the contrast with 
respect to the FLRW spacetime to which the LTB model converges asymptotically. 
See figures 4, 6b, 8b, 10b and lib. The exceptional case is with spatial curvature 
perturbations when this FLRW spacetime is Einstein de Sitter, in which case the 
non-local perturbation 5^ cannot be defined, see section 9). 

We examine below how these differences relate to the shape of the radial profile of 
a scalar at an arbitrary 3 T[t] and the signs and magnitudes ("amplitude") of the 
perturbations. 



7.1. Signs of the perturbations vs type of radial profiles. 

Consider a scalar A whose radial profile is monotonic throughout a given domain [^] 
Then, following 1(23), the clump/void profiles are defined by 

Clump profile A' < =>> A' q < => A < A qi (34a) 
Void profile A' > => A' q > A > A q , (346) 

where A(0) = A q (0) = A c and A'(0) = A' q (0) — and we have assumed absence of 
shell crossings and singular layers (hence R' > holds everywhere or it has a common 
same order zero with A' and A' q |T51 [H]). The relation between the signs of the 
perturbations and the profile type follows from ^ (it can also be seen schematically 
in figures 1-4 and in the graphs of the profiles and perturbations of figures 6-11 of 
the numerical example of section 6): 

• Local perturbations: 

! clump profile : 5^ < 0, 
(35) 
void profile : 5^ > o. 

* Such domains always exist for n finite, even if A is not monotonic for values r > r^. The nature 
and evolution of radial profiles of LTB scalars were examined extensively in reference 1191 . In general, 
density void profiles are fully compatible with regularity conditions (absence of shell crossings) for 
hyperbolic models, but not for elliptic or parabolic ones. Also, for models with a non-simultaneous 
big bang (t' bb ^ 0, nonzero decaying modes) density void profiles always emerge from initial clump 
profiles after a transition ( "profile inversion" , sec figures 6a, 6b and 7a) , whereas a density void profile 
exists for all the time evolution only in models with a simultaneous bang (zero decaying modes). A 
profile inversion of the expansion scalar H necessarily occurs in elliptic models but not in hyperbolic 
ones (see 1191 for further detail). 
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Figure 4 depicts this sign relation. It can also be appreciated by comparing the 
profiles of p and T-L in figures 7 and 10 with the corresponding local perturbations 
(JW and 6™ in fi gures 6a, 8a, 9a and 11a. 
Non local perturbations: 

Finite domains T>[r b ] and Swiss cheese holes. The sign relations in (34a|) and 



(346 1 imply that there always exists a value r — y < r b such that A{y) = A q (r b ) = 
(A} q [r b ], and thus we have from (12): <5nl (y, r b) = 0. Since A c > (A) q [r b ] holds 



in clump profiles and A c < (A) q [r b ] in void profiles, then 6^ changes sign as 
follows (see panels (b) and (d) of figures 2 and 3): 

Clump: 5$ > for < f < y, 6$ < for y < f < r b , (36a) 

Void: Sif < for < f < y, 4l } > for y < f < r b , (366) 

Asymptotic background In the limit r b — > oo we have A — > A^ , (A) q [r b ] — > A^ 
and T>[rt,] — > 3 T[t], hence for all r: 

Clump: A > A x 5$ > 0, (37a) 

Void: A < A x 6$ < 0, (376) 

The relation between the type of profile and the sign of the perturbations is 
illustrated schematically by figure 4. It also emerges by comparing the profiles 
of p and H in figures 4, 7 and 10 with the signs of the local and non-local 
perturbations an d $^ in figures 6b, 8b, 9b and lib. 



Comparing the sign relations in ( |35| ) with those in (36a)-(366) and (37a)-(376) (and 
looking at figures 6b, 7a and 10b) clearly indicates that non-local density perturbations 
^nl* have the expected intuitive signs that we tend to infer for perturbations through 
the familiar notion of a "density contrast" with respect to a fiducial (or background) 
FLRW density value p: an over-density (clump) is a positive perturbation because 
local density is larger than p, while an under-density (void) is a negative perturbation 
because p is smaller than local density, which we can now identify with (p) q [rb\ or its 
asymptotic limit (see figure 6b). In fact, this intuitive notion of a "contrast" can 
be applied to all covariant scalars, not just the density (it is applied to the Hubble 
scalar H in the numerical example in figures 8b and lib). On the other hand, the 
"gradient" based local perturbations 6' > have the opposite sign to the intuitively 
expected one and thus are more abstract and counter-intuitive (compare the profiles 
of p, % in figures 7 and 10 with the signs of 5^ p \ 5^ in figures 6a, 8a, 9a and 11a). 

7.2. Amplitudes. 

The difference between the perturbations 5^ and #4l can also be understood in 
terms of the notion of "amplitude" often introduced in the context of simple intuitive 
Newtonian [51] and relativistic [6l [HJ [12] perturbations. In order to illustrate this 
point we consider an asymptotic perturbation in an LTB model converging radially 
to FLRW (as in figure 4), so that for all A we have along the slices 3 T[t) the 



asymptotically FLRW behavior of (31 ) 



lim A(r b ) = lim (A) q [r b ] = A^t) = A(t) ^ 0, (38) 



n- 



where A(t) is the A scalar for the asymptotic FLRW state (the exception is the spatial 
curvature if the FLRW state is Einstein de Sitter, as in this case JC^ = K = 0). 
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Figure 3. Confined perturbations. The panels display the local and non- 
local perturbations (<5^' and ) f° r * ne radial profiles displayed in figures 1 
and 2. Notice that both types of perturbations vanish at r = ri, when there is a 
Swiss cheese matching (panels (b) and (d)). Notice also that <5'^' and have 
opposite signs around the center. 



Assuming for simplicity a monotonic profile, then without loosing generality we can 
express any A in any complete slice 3 T[t] (not intersecting a singularity) in terms of 
a "contrast amplitude" (f>A as 



A = A x [1 + A f(r)} , with ^ = 4f(0,oo) 



lim 



A r 



(A) q [n] 



Ac 
A 



(39) 

where A c = A(t, 0) = A q (t,0) and / is a smooth non-negative function satisfying 
fc = /(0) = 1 and /'(0) = 0, as well as f(r), f(r) — > as r — > oo. Since we can 
always choose at any 3 T[t] the radial coordinate as R = r, so that R 2 R'dr = r 2 dr, 
the form (39) implies from 1(13) and 1(18) 

J r /f 2 df 



A q (r)=A oo [l + cj> A f g (r)}, f q (r) = (f) q [r 
while ([8]), (pi and p4| yield 

/ 5M = 



Jo r2df 



(/ - u 



(40) 



(41) 



l + 4>Afq ' 

which illustrate the following basic features: 

• The contrast amplitude <$> A obeys a simple linear relation with the more intuitive 
non-local perturbations S^t\ but its relation with local perturbations 5^ is 
non-linear. This is consistent with the fact that the latter perturbations are less 
intuitive than the former. 
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Figure 4. Asymptotic fluctuations and perturbations. Panels (a) and (b) 
correspond to a clump profile, while ( c) and (d) depict a void profile. Since the 
domain encompasses a whole time slice 3 T[t], non-local fluctuations compare local 
values of A(r) with the global average (A) q [?•;,] — > A x in the limit — > oo, which 

is equal to the corresponding scalar A(t) of the FLRW asymptotic state. Non— 
(A) 

local perturbations <5^ L yield the "contrast" between local values of A and this 
average. The amplitude of the perturbation is <pA, while ip^ depicts the maximal 
value of the local perturbation 5( A K Notice that 5^., ' is positive/negative for the 
clump/void, while S^ A ^ has opposite signs. Also, both perturbations vanish as 
r — > oo. The graphs of asymptotic perturbations displayed in figures 6, 8, 9 and 
11 corroborate the qualitative forms shown in panels (b) and (d). 



• Both perturbations vanish in the domain boundary r —¥ oo where the 
correspondence with the FLRW background occurs, but their behavior at the 
center is different: 6^(0,00) = <f> A ^ (in general) while 5^(0) = (see figures 
4, 6, 8, 9 and 11). 

• The relation between the signs of the perturbations and the type of profile can 
also be given in terms of the amplitude: 

Clump profile: A c >A x , 4> A > 0, 8$ > but 6 {A) < 0, (42a) 

Void profile: A c < A x , <\> A < 0, 8$ < but 8 iA) > 0. (426) 

wher e we u sed the fact that (by construction) /' < so that the sign relations 
( 34a H 346) imply / — f q < and the sign of S^ A ' is the opposite of the sign of 
4>A- 

Evidently, the convention whereby a clump or void {i.e. overdensity or underdensity 
when A — p) respectively correspond to positive or negative contrast follows from the 
sign of the amplitude contrast embodied in $nl\ cind is more intuitive and easier to 
follow than the complicated relation between radial gradients A 1 associated with 5^ K 




Figure 5. Amplitudes of local and non-local perturbations. The graph 
depicts the monotonic non-linear relation between the amplitude, ipA , of non-local 
perturbations and the maximal value of the local perturbation, tpA ( see figure 
4), which shows that both perturbations provide a self-consistent qualitatively 
analogous description of the deviation from FLRW. Only for small amplitudes 
\4'a\ ^ 1 the relation is linear. 



However, in spite of their differences in signs when referred to clump/ void 
radial profiles, the deviation form homogeneity (FLRW conditions) can be traced 
effectively with both types of perturbations, S^} and 5^ A \ because there is a consistent 
monotonic relation between their magnitudes (their maximal/minimal value in a given 
domain) and the amplitude 4>a ■ For domains with monotonic / the the amplitude cf> A 
of the non-local perturbations 5^} coincides with their extremal value, which occurs 
at the center r = (see figures 6b, 8b, 9b and lib). For local perturbations the 
relation between 4>A and the extrema of 5^ (denoted by ip A , see figure 4) is more 
complicated, as the latter may occur for different values of r for different slices 3 T[t] 
and depends on the choice of / (see the maximum/minimum of the local perturbations 
in figures 6a, 8a, 9a and 11a). Since a result for a general / is hard to find, we consider 
the special form / = (1 + r™) _1 where n > (the results are qualitatively analogous 
for all other forms). The result, shown in figure 5, illustrates the non-linear monotonic 
relation between ip A and 4>A- Hence, both types of perturbations provide a consistent 
estimate of the deviation from homogeneity of the models in which larger values of 
\(/)a\ correspond to larg er magnitudes of and 

8. Linear limit. 

While the 5^ and the are exact quantities that need not be "small' and do not 
comply with linear evolution equations, we can expect (intuitively) that they should 
somehow reduce to linear dust perturbations (in the comoving gauge) when their 
magnitudes (i.e. amplitudes) are small, as in the early times regime in the example 
of section 9. A more rigorous way to examine their connection to linear perturbations 
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of dust sources follows by constructing second order equations for and 5^'- 



• Local per turb a tion s. Wc differentiate both sides of ( 19a) and use the remaining 



equations (19&)-(19d) to eliminate all derivatives except Hr- p > and S^ p \ leading to 
2 [j^] 2 



1 + 5(p) 

• Non local pert urbations. 



Considering asymptotic backgrounds (rf, 
2 [6 



+ 2H q <5 W - 4irp q <yW (l + 5^ = 0, (43) 

We follow the same procedure as above: 
ewrite it in 
oo ) yields: 



differentiate &25ah and use (j25&h — d26 £>h to rewrite it in factors of J^l and <*>nl 



■'NL 



1 



dp) 



•4(67^ -AH q )S } 



(p) 

NL 



6 (H q - H, 



1 + & 



= 0, 
(44) 



which holds for confined domains by replacing p^, with (p)q[rb\, CH)q[fb] 



We remark that both (43) and (44) are exact non-linear equations for #W or <5 



(p) 



(equations similar to (43) have been obtained in [521 153] for Szekeres models). Also, 
notice that (43) and (44 1 coincide as r — > oo for asymptotic perturbations and as 
r — > r& for confined domains 2?[rt,]. 

In general there is no reason to assume that the amplitudes (\>a or their time 
derivatives 4>a are small in a generic LTB model subjected to non-linear evolution 
equations, but if in a given range of t we have \4>a\ ^ 1 then the relation between 
both 6^ and 5^} with 4>a should become similar. This can be seen by expanding for 
\4>a\ <^ 1 in the non-linear relation in (41) between the 5^ and the amplitude 4>a '■ 

<5 (A) «(/-/ 9 )0a + O(^), (45) 

so that the become linear on <j>A (as the 8^}), and thus the magnitudes of 8^ 
and 8^ are both proportional to <Pa because < / < 1 and the fluctuations f — f q are 
bounded. As a consequence, if \4>a\ "C 1 both perturbations are of the same order in 
4> A1 and thus both comply with \8^\ <C 1 and \8^\ <C 1 and the time evolution 
of 5^ must be the same as that of Si^} at leading order in <p A . In particular, 
if we take A = p, H, the linearity conditions \<f> | <C 1 and \<j) H \ <C 1 imply that 

\5V>\ « 1, \6$\ « 1 and \6<M\ « 1, |^| « 1 hold, and thus « 1, |j<#| « 1 
also holds (from (19c) and ( |25a[ )). Therefore, we also have p q rj p^ rj p and 
Hq « ~ % (the same relations hold for finite domains by replacing , 



with {p} q [ro}, (H) q [ro\). Considering all these implications, if \<Pa\ *C 1 the second 
order evolution equations ( 43 ) and ( 44 ) becomes at order 4> A 

S(p) + <5 (p) - 4tt Px = 0, (46) 

which is formally identical to the evolution equation for gauge invariant linear density 
perturbations of a dust source around a FLRW background characterized by p^ , 
(or (p) q [rb\, (H) q [rb\ for bounded domains) in the comoving gauge [5TJ[5I], which for 
dust is a synchronous gauge as well. 



9. A numerical example of a cosmic density void. 



In order to illustrate the utility of the evolution equations of section 4 for model 
building, as well as the properties and differences between local and non-local 



Weighed scalar averaging in LTB dust models, part II 



22 



perturbations depicted qualitatively by figures 4b and 4d, we consider the evolution of 
a cosmic void configuration that emerges from the numerical solution of the systems 
fll9ol)-([l9^1 and ([25a|)-p6b| [| . For this ] purpose, we consider a regular hyperbolic 
model (JC q < or < fl q < 1) that is radially asymptotic to a spatially flat (JC = 0) 
FLRW model (Einstein de Sitter) [THj characterized by the Hubble factor, density and 
big bang time H, p, t bb that satisfy 8np/(3II 2 ) = 1 and H(t — t bb ) — 2/3 for all t. 
Hence, initial conditions must be selected such that shell crossings do not arise (see 
Appendix C of [TI5]) and the limits K. q — > or ft q — > 1 as r — > oo hold along every 
time slice 3 T[t] [H]. 

We consider the last scattering surface as the initial time slice, hence the subindex 
o will correspond to evaluation at t = to = £ LS , which suggest using the constant 
H~<} = if _1 (t LS ) as the characteristic length scale. The nearly homogenous and 
spatially flat conditions at £ LS imply that initial value functions must satisfy: 



H 



qO 



n 



qO 



H^S, PqO 



3H q0 



1. 



%go(*LS 



tu 



tc 



qO 



0. 



~ 2 

^bb] ~ g 1 



(47a) 
(476) 



for all r, together with the strict limits p q o — > /3 LS , H q Q — > Hl S , ^qo —> 1 and t bb — > t bb 
as r — > oo. By normalizing with respect to i? LS we obtain the following dimensionless 



initial value functions that are compatible with (|47a|)-(|476|) 



f-qO 



l q0 



HqO 

Ht.f, 



1 



1 + (r/n) n 



K q0 



c 



qO 



'^qOI 



il/2 



1 + (r/r 2 )P 

-,1/2 



(48 a) 



(486) 



l + (ryYi) n l + (r/r 2 )P_ 

where the positive constants e, Cj n, p comply with e <C 1, C *^ 1 } n < 3, p < 2 
(the latter two follow from regularity conditions |18j). while r%, r 2 are constant length 



scales, and we used the constraint (66) 



Considering the initial value functions p q o, h q o in (48a)-(486) for the parameter 
values n — 3, p = 3/2 , r\ = 1/2, r 2 = I/a/3, e = 0.001, ( = 0.01, we integrate 
the system (19a)-(19d) for the dimensionless variables p q = 4irp q /(3Hl s ), h q = 
Hq/H^s, 8( p >, 8^ in terms of the dimensionless time 

t = H LS (t - t bb ) - |, (49) 



so that t = corresponds to t = i LS , since t bb ~ t bb with t bb — t bb (0) ~ 0(max(e,£)). 
The bang time is marked by r « —2/3 + 0(max(e,£)), though the LTB model is 
no longer valid for t < t LS (or —2/3 < r < 0), and thus we only consider the 
range r > 0. Since if LS = H (l + z) 3 / 2 and z(t LS ) ~ 1100, the present value of 
Ho 70km/(sec Mpc) yields for the present time the value r = 3.6 x 10 4 , which 
approximately corresponds to 13 Gys. 

By considering non-local perturbations that are asymptotic (i.e. rj, — > oo taking 
the form p3|), t ogether with a FLRW background state that is Einstein de Sitter, 
the system ( 25a)— ( 266 1 simplifies considerably, since the terms (p) q [rb\, (H) q [rb\ in 



these equations take the asymptotic forms p, H that comply with (Q) Cj 



1. 



tl We do not claim that this void model is "realistic" nor that it provides a good fitting to observations. 
Its purpose is simply to illustrate the evolution of the local and non-local perturbations of p and H 
and its relation to the radial profiles of these scalars. 
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(a) 



(b) 



Figure 6. Radial profile of density perturbations (early times). Panels 
(a) and (b) respectively depict local and non— local perturbations for values 
of r corresponding to early cosmic times rj, i = 1,2,3,4 corresponding to 
0, 0.5, 1.0, 5.0. The profiles have the same form as depicted in the qualitative 
diagrams of figures 4b and 4d. The thick curves mark the last scattering surface 
t — 7~i — 0. Notice how the perturbations invert their sign between n and T2, 
corresponding to the density profile inversion shown in figure 7a (see the text for 
further explanation). 



Hence, all terms involving (p) q [rb\, (H) q [rb\ can be replaced by closed exact analytic 
dimensionless forms p,, h associated with the Einstein de Sitter asymptotic state, which 



can be given analytically in terms of the dimensionless time (49 1 by: 

~ H 1 47T/5 h 

ZHh 2 



1 



1 

253' 



-i 2/3 



1 



(50) 



However, from a computational point of view, it is easier to obtain the non-local 



perturbations through the relation ( 24 ) 



fl \ 



1 + tfOO 



1. 



(51) 



where jl, h are given by (50) and [i qi h q , S^ p \ 5^ follow from the numerical 
integration of (19a)-(19rf|. 



9.1. Early cosmic times. 

Since At — 1 corresponds to a time lapse H^, 1 ~ 10 5 ys, values r ~ O(l) 
should describe a linear regime with small amplitudes for both local and non-local 
perturbations. This can be appreciated in figures 6 and 8, depicting the radial profiles 
of local and non-local perturbations of the density and the Hubble scalar for the values 
[n, r 2 , r 3 , r 4 ] = [0, 0.5, 1.0, 5.0] that range between the last scattering surface (r = 0, 
marked by thick curves) to about 10 6 ys (r = 5). 

The difference between local and non-local perturbations and its relation to the 
perturbation amplitude and type of radial profile (clump or void) clearly emerges by 
comparing density perturbations in figures 6a and 6b and the radial density profiles 
in figure 7a. Notice how the amplitudes of the perturbations, (f) p and ip p , exhibit the 
behavior described in section 7 and depicted in figure 4. As shown in figure 7a, the 
initial configuration is a clump with small amplitude (thick curve marked by r = 0) : it 
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0-5 1 1-5 2 2.5 o 0.5 I 1.5 2 _ 

(a) (b) 

Figure 7. Radial profiles of the density and Hubble scalar (early times). 

Radial profiles of the local scalars p and "H. normalized by their central values 
p c = p(t, 0) and W c = W(t,0) for the same values of r of figure 6. The profiles 
have the form sketched in figures 4a and 4c. Notice the transition from a clump 
to a void profile of p between ti and T2, while "H has a clump profile for all T. 




0.5 I 1.5 2 f 0.5 1 1.5 2 f 

(a) (b) 

Figure 8. Radial profiles of perturbations of the Hubble scalar (early 
times). Panels (a) and (b) respectively depict local and non— local perturbations 
for the same values of r of figure 6. The sign of the latter is characteristic of a 
clump profile (see figure 7b). 



is a positive non-local initial perturbation (i.e. positive density contrast) depicted by 
the thick curve in figure 6b, while the initial local perturbation (thick curve in figure 
6a) is negative (negative radial gradient of the density). As the evolution proceeds, 
both types of density perturbation in figures 6a and 6b have reversed their sign at 
t = r 2 , with the profile passing from that of a clump to a void in figure 7a [ff] : 
the non-local perturbation becomes negative (voids have negative density contrast) 
and the local one becomes positive (radial gradient of density is positive in voids). 
However, the perturbations of the Hubble scalar in figures 8a and 8b do not change 
sign: the non-local perturbation is positive (positive contrast) and the local one is 
negative (negative radial gradient), and thus the radial profile is that of a clump for 
all these values of r (see figure 7b). 

ff The occurrence of this profile inversion in regular hyperbolic models is reported in |19) . A necessary 
condition for it is 6$ - (3/2)5 ( q ) > 0. 
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2 .1 4 f, 5 0.5 1 1 .5 2 f 



(a) (b) 

Figure 9. Radial profile of density perturbations (late times). The 

panels respectively depict local and non— local perturbations for time values 
i = 1,2,3,4,5,6, 7,8,9, 10 given by 100, 1000, 5000, 18000, 36000, 72000, with 
the thick curve marking present day cosmic time Tg = 36000, which roughly 
corresponds to 13 Gys. Local perturbtions have positive sign while the non-local 
ones are negative, all of which is characteristic of void profiles (see figure 10a). 
The amplitudes clearly correspond to the non— linear regime. 




Figure 10. Radial profiles of the density and Hubble scalar (late times). 

The radial profiles of the same scalars of figure 7 for same cosmic times as figure 
9. Notice the density void profile, while the Hubble scalar keeps a clump profile. 



9.2. Late cosmic times. 

We consider now cosmic times 75, t 6 , t-j, t 8 , 79, r 10 respectively corresponding to 
[100, 1000, 5000, 18000, 36000, 72000]. Hence the curve marked by r 9 = 36000 
approximately corresponds to present day cosmic time and is depicted by thick curves. 
The radial profiles of local and non-local density perturbations are displayed for these 
values of r by figures 9a and 9b. Evidently, the perturbations' amplitude increases 
to large values that indicate a non-linear regime and the density profile is that of a 
density void (figure 10a). The sign of the non-local perturbations is negative, as the 
density contrast is negative in voids, reaching an over ~ 90% negative contrast at the 
center for late cosmic times (as in the void models of chapter 4 of [12]). However, local 
perturbations are positive, as radial gradients are positive in void profiles. On the other 
hand, as shown by figures 11a and lib, the local and non-local perturbations of the 
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Figure 12. Radial profile of Q q and spatial curvature perturbation S^K 

Panel (a) displays the radial profile of the q— scalar fl q , defined by j7|. The 
curves shows the transition from nearly homogeneous spatially flat conditions at 
t = t L s (Q q is 1) to i very inhomogeneous void with negative curvature with 
SI, <C 1 at the void center in the present day cosmic time. Panel (b) depicts the 
local perturbation of the spatial curvature 8^'. Since K. — > as r — > oo, this 
perturbation is large with 8^ ' — > —1/2 in this limit. 



Hubble scalar are, respectively, negative and positive, characteristic of a clump profile 
(figure 10b). While the amplitudes of the local perturbations of the Hubble scalar 
remain small and within linear regime (~ 1CP 2 ), the amplitudes of these perturbations 
show a steep growth from figure 9a to 11a. On the other hand, the amplitudes of the 
non-local perturbations of the Hubble scalar become large, showing a present day 
~ 30 % positive contrast. 

The radial profiles of the q-scalar Vt ql depicted in figure 12a, describe the 
evolution of this LTB model as a transition from nearly homogeneous and spatially 
flat conditions (Q q ~ 1) at t = £ LS to a very inhomogeneous void with Q q <C 1 with 
marked negative spatial curvature. Since the negative spatial curvature complies with 
IC q — > for all t, the non-local spatial curvature perturbations cannot be defined 
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because {K) q — > K. = 0, but there is no problem in defining local perturbations 5^'. 
As shown in [T8] , the local perturbat ion o f a sca lar tending asymptotically to zero is 



not small. For the initial conditions ( 48a H 486 1, $^ reaches an asymptotic limiting 
value -1/2 as r — > oo (see figure 12b). 

10. Conclusion and summary. 

We have examined the perturbations that emerge from the weighed averaging 
formalism developed in part I. These perturbations can be local (5^ A ') or non-local 
(4l j! depending on whether they are, respectively, defined for q-scalars that are 
functions or functionals (a summary of part I is given in section 2). By introducing 
an initial value parametrization, we showed in section 3 that the q-scalars common 
with FLRW models (A g and (A) q [r b ] for A = p, H, JC, fi) identically satisfy FLRW 
scaling laws that mimic FLRW expressions commonly used in LTB void models that 
probe the possibility of explaining observations without resorting to dark energy 
[Ml H3 113 HH H3 1301 EQ El [MIIMIISS]. Also, the q-scalars determine (via Darmois 
matching conditions) a unique FLRW background state that can be either a formal 
reference FLRW model that is different for each domain, or an actual background 
FLRW spacetime if a domain 2?[r&] is smoothly matched to a FLRW region (Swiss 
cheese models). The following results are worth highlighting: 

• As proven in part I (see I(42a)-I(42b) and 1(43)), the local density perturbation 
5( ps> is directly expressible in terms of the ratio of Weyl to Ricci curvature: ^2/^, 
where \&2 is the only nonzero Newman-Penrose conformal invariant and TZ is the 
Ricci scalar, whereas 5^ is expressible in terms of the ratio of anisotropic vs. 
isotropic expansion: £/%, where £ is the eigenvalue of the shear tensor. All q- 
scalars can be given in terms of these scalar invariants by means of the constraints 



(10 Mil I 



The evolution equations for the q-scalars and their perturbations (local and non- 
local) completely determine the dynamics of the models, and thus provide an 
alternative to the use of analytic solutions in the study and applications of the 
models. Although the resulting systems involve PDE's, they contain only time 
derivatives and the constraints are purely algebraic in nature and preserved in 
time. Hence, these equations can be treated effectively as ODE's (section 4). 
We tested the numerical integration of these systems in the example provided 
in section 9 of a void model that is radially asymptotic to an Einstein de Sitter 
FLRW model. 

Since the q-scalars behave effectively as FLRW scalars and their perturbations 
convey for every domain the (exact) deviation from FLRW dynamics, we can 
rigorously re-interpret LTB dynamics as the dynamics of exact spherical dust 
perturbations on a FLRW background defined by the q-scalars. As mentioned 
before, this background state is an abstract reference background defined at every 
domain by the continuity of the q-scalars under Darmois matching conditions 



(16a|-(16c). It is an actual FLRW spacetime only if a smooth matching with a 
FLRW region is considered in the context of "Swiss cheese" models (section 5). 

Both local and non-local fluctuations and perturbations can be either confined in 
a given domain I>[r&] (with or without assuming a "Swiss cheese" configuration 
through a matching with a FLRW region, see figures 1, 2 and 3), or an asymptotic 
perturbation for the case when T>[ri,] becomes the whole time slice 3 T[t] in the 
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limit rb — > oo (see figure 4). The choice between confined and asymptotic 
perturbations depends on the boundary conditions of the specific problem or 
application that we may work out with LTB models: 

— Confined perturbations (local or non-local) without a matching with FLRW 
are practical because they are easily applicable to study the inhomogeneity 
pattern of any generic LTB configuration, but are more abstract because the 
FLRW background state is a fictitious reference spacetime that changes for 
each domain. 

— Swiss cheese configurations are more intuitive because the background state 
is a single (and non-fictitious) FLRW spacetime. They are useful if we wish 
to describe various inhomogeneous regions, as in Swiss cheese void models 
used to fit observations J5B S 29, 30, 3 1] (see reviews in and [35]), but their 
intuitiveness is offset by the emergence of artificial jointly packed layers near 
the hole boundary ("humps" and "bags" in the radial profile (see figures 1, 2 
and 3) that arise because of need to impose continuity of the perturbations, 
which must vanish in a non-fictitiuos background located at a finite radius 
(also, shell crossings necessarily arise if the LTB region is hyperbolic, see 
Appendix A). 

— Asymptotic perturbations do not present these inconveniences, and thus are 
more natural, and also practical if the radial profile of the scalars rapidly 
converges to a given set of FLRW values (as is the case in the model presented 
in section 9 and illustrated by figures 6-12). 

The local perturbations provide a measure of inhomogeneity through the 
local magnitude of the radial gradients of covariant scalars, while non-local 
perturbations do so through the familiar notion of the "contrast" between local 
values of these scalars and a reference FLRW value defined as a q-average for 
each fixed whole domain, or for the whole slice in the asymptotic limit. Even 
if their evolution equations are more complicated, non-local perturbations are 
more intuitive than local ones, as they allow us to associate over and under 
densities with (respectively) positive or negative density amplitudes, while local 
perturbations yield the opposite sign. We have illustrated these points through 
the notion of the "amplitude" of the perturbations (see sections 6 and 7 and figure 
5 and compare the profiles of p, H in figures 7 and 10 with the profiles of their 
local perturbations in figures 6a, 8a, 9a and 11a). As a possible application, the 
exact local and non-local perturbations may be useful in understanding structure 
formation through an approach that involves the radial gradient and profile of the 
perturbations (as for example in [55]). 

The definition of non-local perturbations 6^ by the maps ( 29 1 and ( 30 1 provide 
(through the q-average and its relation to kinematic and curvature invariants) 
a rigorous and coordinate independent interpretation for simple examples of 
perturbations that are based on the notion of a contrast with respect to a FLRW 
background. These simple contrast perturbations are frequently introduced as 
ansatzes in many text-book and articles, for example, in the context of simple 
Newtonian models of structure formation (the "top hat" or "spherical collapse 
model" [51]) but also in linear perturbations [54] and in astrophysical applications 
of LTB models (see examples in [6j [HI [12] ) . 

Both local and non-local perturbations yield in the linear limit the familiar dust 
perturbations of linear theory in the comoving (or isochronous) gauge (section 8). 
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• The example of a void model given in section 9 (which follows from the numerical 
integration of the evolution equations of section 4) clearly illustrate the properties 
of local and non-local perturbations described above. It also shows the potential 
utility of the q-scalars and their perturbations as tools for LTB model building. 

While the q-average formalism cannot be used to study and understand the 
dynamical relation between perturbations and back-reaction, it does provide through 
the results of part I and the present article (part II) interesting theoretical connections 
between averaging, perturbation theory, invariant scalars and statistical correlations 
of p and H, which signals a valuable theoretical insight on how the averaging process 
should work in any generic solution of Einstein's equations (at least in LRS spacetimes 
whose dynamics is reducible to scalar modes). A coordinate independent study and 
generalization of the "growing" and "decaying" models is still necessary for a 
complete theoretical study of the perturbations that we have examined. This task 
will be undertaken in a separate article currently under preparation. 

Since most formal and theoretical results obtained for LTB models can be readily 
applied to Szekeres models |23| . the same type of exact perturbation formalism 
can be devised for these non-spherical models and perhaps even to more general 
spacetimes. Specifically, we propose identifying fluid flow scalars that comply with 
FLRW dynamics (hopefully related to some weighed average), and then expressing 
local 1+3 scalars as fluctuations of the new scalars. The following step would be to 
explore if the 1+3 evolution equations given in terms of the new set of variables has 
the structure of evolution equations and constraints in the context of a formalism 
of exact perturbations on a FLRW background. We feel that this proposal is worth 
considering in future research. 
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Appendix A. Shell crossings and Swiss cheese models. 



In a Swiss cheese model in which the FLRW background is a spatially flat Einstein 



de Sitter model the second of the Darmois matching conditions (16a|-(166) and the 



supplementary condition (17) must be replaced by: 

qb = (K) q [fb] = ic(t) = o, n qb = (n) g [r b ] = o(t) = 1, s, 



K 



(X) 



0. 



(A.l) 



where r b is bounded. We prove in this Appendix that these conditions are incompatible 
with the following condition to prevent shell crossings in hyperbolic models JC q o < 
(see [7] and Appendix C of [H?]): 

rE' 2 2 



3E 



2 

> 
3 ~ 



3 ~ 



E' = (|/C 90 |r 2 )' > 



(A.2) 



In order to comply with (A.l | we must have \IC q o\ — > as r — > r b , and thus |/C 9 o|' — > 



must also hold in this limit, but (A.2) requires l/Cqor to be monotonically increasing 
in the range < r < r b and we have E = \K q0 \r 2 -)• as r -> [?1 IT U H5 1 ITS). 
Therefore, E' necessarily reverses its sign in this range. As a consequence, shell 
crossings necessarily occur in all Swiss cheese models with an Einstein de Sitter 
background in which the LTB section is hyperbolic (regardless of the type of profile 
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of the density). This problem, which has been reported by [26] and commented in 
chapter 5.3.5 of [12], does not arise if the LTB region is elliptic (as the sign of E' is 
not constrained) and in hyperbolic models radially converging asymptotically to an 
Einstein de Sitter state (as the void model of section 9), since in this case functions 
JC q o can be found such that both \fC q o\ — > as r — > oo and (\K. q o\r 2 )' > hold [18 : 
for example, if \JC q o\ decays as r~ p with < p < 2 (as we assumed in (48a|). 
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